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ABSTRACT 

Oceans general circulation models predict that global warming may cause a decrease in the 
oceanic 0 2 inventory and an associated 0 2 outgassing. An independent argument is presented 
here in support of this prediction based on observational evidence of the ocean’s biogeochemical 
response to natural warming. On time scales from seasonal to centennial, natural 0 2 flux/heat 
flux ratios are shown to occur in a range of 2 to 10 nmol 0 2 per Joule of warming, with larger 
ratios typically occurring at higher latitudes and over longer time scales. The ratios are several 
times larger than would be expected solely from the effect of heating on the 0 2 solubility, 
indicating that most of the 0 2 exchange is biologically mediated through links between heating 
and stratification. The change in oceanic 0 2 inventory through the 1990's is estimated to be 
0.3X0.4P10 14 mol 0 2 per year based on scaling the observed anomalous long-term ocean 
warming by natural 0 2 flux/heating ratios and allowing for uncertainty due to decadal variability. 
Implications are discussed for carbon budgets based on observed changes in atmospheric 0 2 /N 2 
ratio and based on observed changes in ocean dissolved inorganic carbon. 

INTRODUCTION 

Repeated hydrographic surveys indicate that the upper 3 km of the oceans have warmed 
(1) and intermediate waters of high-latitude origin have freshened (2) over the past few decades. 
Model studies indicate that upper-ocean warming, high-latitude freshening, and an associated 
increase in the density stratification of the upper ocean are expected consequences of the changes 
in atmospheric radiative forcing caused by fossil-fuel burning and other human activities (3,4,5). 

Repeated hydrographic surveys also indicate that small detectible changes have occurred 
in oceanic dissolved 0 2 concentrations. As summarized in Table 1, detectible decreases in 0 2 
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have been found in intermediate waters in the North Pacific, North Atlantic, South Pacific, and 
South Indian oceans, while small increases have possibly been found in deeper waters in the 
North Pacific and South Indian Oceans. What caused these O 2 changes is unclear, and different 
mechanisms, including changes in ocean circulation rates (6, 7, 8, 9), changes in preformed values 
(10), and changing Redfield ratios (11), and changes in biological production (8) have been 
offered as possible explanations in different regions. While the changes may partly reflect 
natural decadal variability, the clearest 0 2 changes, found at intermediate depths, are in the 
direction of decreasing 0 2 concentrations. A global reduction in dissolved 0 2 is predicted by 
ocean general circulation models (OGCMs) driven by increasing greenhouse gases 
(12,13,14,15). In the model simulations, most of the 0 2 decrease is attributed to enhanced 

stratification. 

Stratification has two competing effects on subsurface oxygen concentrations. First, it 
reduces the transport of nutrients from deeper waters into the upper ocean, thus decreasing 
photosynthetic production and the associated flux of organic detritus into the ocean interior. 
This flux is often referred to as the “biological pump” (16), and reducing the rate of this pump 
increases subsurface 0 2 concentrations by reducing subsurface 0 2 utilization rates. Second, 
stratification limits the downward transport from 0 2 from well-oxygenated surface waters into 
the ocean interior, which serves to reduce subsurface 0 2 concentrations. 

In the modeling studies cited above, the effect of stratification on 0 2 transport exceeds 
the effect on subsurface 0 2 utilization, leading to a net 0 2 decrease. This result is expected 
considering that stratification allows for more complete biological utilization of nutrients in 
surface waters, thus lowering the “preformed” (i.e. initial) nitrate and phosphate content of 
waters sinking into the oceans interior. To conserve the total ocean nutrient inventory, an 
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increase must occur in the inventory of non-preformed nutrients, i.e. those which accumulate in 
subsurface waters from oxidative decomposition of organic detritus. Since 0 2 is consumed by 
organic decomposition in proportion to the amount of nitrate or phosphate released (17), 
subsurface 0 2 inventories can be expected to decrease in response to increased stratification. In 
general, the competing effects of the biological pump and vertical mixing on 0 2 concentrations 
can be assessed based on their net impact on the vertical nutrient distributions. This net impact 
we refer to as the “efficiency” (as opposed to “rate”) of the biological pump. 

Stratification can also be expected to induce a net release of 0 2 from the ocean to the 
atmosphere. The oceanic 0 2 inventory (I, moles) and sea-to-air 0 2 flux (Z, moles/yr) are linked 

according to 

dl/dt = -Z + X dCorg/dt 0) 

where C or g is the oceanic inventory of organic carbon, including both particulate and dissolved 
forms, and X is the 0 2 :C oxidative ratio for destruction or production of marine organic matter, 
and where small terms related to river and sediment transports have been neglected. The term 
X dCorg/dt effectively accounts for the column-integrated net production of 0 2 by marine 
photosynthesis and respiration. This term is presumably small, although not necessarily 
negligible (18), because the main effect of marine photosynthesis and respiration on time scales 
of years to centuries is to redistribute inorganic materials within the ocean, rather than to cause 
accumulation or destruction of the organic carbon (see also ref. 14). If stratification reduces 
subsurface 0 2 inventories via changes in the efficiency of the biological pump, a roughly equal 
outgassing of 0 2 from the ocean to the atmosphere is therefore expected. 

Estimates based on OGCMs for oceanic 0 2 outgassing due to global warming lie in the 
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range of 0.2 to 0.7H0 14 mol yr for the past few decades, with predictions of 1.0 to 1.6 HO 14 
mol yr' 1 for late 21 st century (12,13,14,1 5). According to these estimates, roughly a quarter of 
the predicted outgassing is attributable to the direct effect of warming on the 0 2 solubility, while 
the remainder is due to increased stratification. The predicted changes for the past few decades 
amount to an average decrease in oceanic 0 2 concentrations over a 20-year time frame of 
between 1 and 5 Bmol kg^ 1 if the changes are confined to the top 1000m of the oceans. Larger 
changes are predicted for the Southern Oceans due to reductions in deep convection (12,13). It is 
hard to test these predictions with existing hydrographic data given the sparse coverage and the 
lack of comprehensive syntheses. 

There are several reasons why changes in the oceanic 0 2 inventory could be important. 
Dissolved 0 2 concentration is a useful diagnostic of ocean circulation and biological activity 
which can provide constraints on models of physical and biogeochemical response to climate 
change (13). Small changes in 0 2 content could influence extent of hypoxic regions in coastal 
seas, in sediments, or in the open ocean, with consequences for the cycling of nitrogen and other 
redox-sensitive elements and for the distribution of many marine organisms (19). Changes in 0 2 
are diagnostic of changes in the efficiency of the marine biological pump, which may influence 
the rate at which the oceans absorb anthropogenic C0 2 (12). Finally, at least two approaches for 
estimating sinks of anthropogenic carbon dioxide require corrections for changes in oceanic 0 2 
inventory. The purpose of this paper is to provide an independent estimate of the plausible 0 2 
inventory changes associated with recent global warming and to discuss the implications for 

global carbon budgeting. 
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NATURAL WARMING AND 0 2 OUTGASSING 

It is easily established that a close association exists between ocean warming and 0 2 
outgassing over a range of space and time scales in the open ocean. At middle and high latitudes 
during the spring and summer, when the upper ocean is heated by the atmosphere, the oceans are 
a source of 0 2 to the atmosphere, while in the fall and winter, when the upper ocean is cooled, 
the oceans are a sink of 0 2 from the atmosphere. These seasonal air-sea 0 2 exchanges are driven 
partly by biological exchanges, linked to seasonal stratification, nutrient supply, and irradiance, 
and partly by effects of heating and cooling on 0 2 solubility ( 20 ). 

We have estimated the ratio of seasonal 0 2 outgassing to seasonal heating from global 
archived measurements of dissolved 0 2 in surface waters, climatological winds, and 
climatological air-sea heat fluxes (21). The 0 2 flux/heating ratio varies between 1.5 nmol 0 2 per 
Joule at lower latitudes to 4 or 5 nmol 0 2 per Joule in the 40“® to 60"® latitude bands, as shown 
in Figure 1 . Consistent with Najjar and Keeling (20), we find ratios that are larger than expected 
from the effect of warming on the 0 2 solubility by factors between 1 .5 to 2.5. We further have 
shown, based on comparisons with atmospheric 0 2 /N 2 data (21), that the component of the 0 2 
flux that correlates with heating dominates large-scale seasonal 0 2 exchange. These results 
indicate that seasonal heating, through its effect on stratification, biological productivity, and 0 2 
solubility, is a major driver of the exchange, and not just coincidentally correlated with the 
exchange. 

Linkages between heat fluxes and oxygen fluxes are also evident on longer time scales, 
as revealed from a plot of the tracer 0 2 * = 0 2 + 175P0 4 versus potential temperature (&), as 
shown in Figure 2. The tracer 0 2 , which is identical to Broecker s tracer PO 4 (25) but 
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expressed in 0 2 rather than P0 4 units, is a measure of the 0 2 gained or lost by a water parcel 
through air-sea gas exchange (26). 0 2 is largely conserved below the sea surface, where 
photosynthesis and respiration produce compensating effects on 0 2 and P0 4 . 0 2 keeps track of 
air-sea 0 2 exchanges driven by both solubility changes and by the processes controlling the 
efficiency of the biological pump. Assuming rapid air-sea equilibration, the solubility 
component is simply given by the 0 2 solubility (0 2 ), so that the remaining biological pump 
component is proportional to the (apparent) preformed phosphate content (P0 4 - (0 2 °-0 2 )/175) 
content of the water (26). 

A prominent feature in Figure 2 is the strong association between 0 2 and £7 in waters 
of the main thermocline waters between 6 and 18*®C. Waters around 18^C consistently have 
lower 0 2 * than waters around 6°®C, indicating that conversion of cooler water to warmer water, 
which occurs mostly at low latitudes, leads to outgassing of 0 2 , while the conversion of warmer 
water back to cooler water, which occurs mostly at higher latitudes, leads to uptake of 0 2 . In the 
Indian and Pacific ocean, the 0 2 /l& slope of ~22 Smole kg^'^C 1 in waters between 6°® and 
18^>C, is equivalent to ~5 nmol of 0 2 per Joule of warming, as derived by multiplying by 
seawater density and dividing by heat capacity. The slope is several times larger than expected 
from solubility changes alone, indicating that the 0 2 exchanges are mainly controlled by 
variations in the efficiency of the marine biological pump. A generally weaker and but more 
variable 0 2 /heat relationship is indicated for waters warmer than about 18^C, which is 
consistent with the lower nutrient content of these waters and a reduction of influence of 
biological relative to solubility effects. 

A less steep O 2 */£7 trend is seen in Figure 2 for the thermocline of the North Atlantic 
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ocean compared to other oceans, possibly owing to the lower nutrient content of the North 
Atlantic. The North Atlantic trend of ~13 Bmole kg'^C' 1 (equivalent to ~3 nmole J' 1 ) 
connects low-latitude surface waters in the Atlantic with North Atlantic Deep Water (NADW), 
which lies below the O 2 V 1 & trend of the other oceans. A steeper O 2 /& ratio of around 30 to 
40 Hmole kg^ l= ®C 1 (equivalent to 7.5 to 10 nmol J ') is found for the deep Antarctic sequence 
(25), which is driven by ventilation of deep waters around Antarctica. Here stratification 
induced by warming and freshening in the summer months inhibits the uptake of 0 2 by deeper 
waters, while the breakdown of stratification induced by winter-time cooling and brine rejection 
from sea ice enhances O 2 uptake and deepwater formation (27). The air-sea exchanges and 
watermass-mixing around Antarctica effectively convert circumpolar deep water, which derives 
largely from NADW, into colder Antarctic surface waters. By this conversion, the oceans 
around Antarctica release heat to the atmosphere and take up O 2 . 

Due to sparse coverage, Figure 2 doesn’t resolve well water masses in the Equatorial 
Pacific. Here, however a very different relationship between heating and 0 2 flux is known to 
exist. Equatorial upwelling raises cool, oxygen-deficient waters to the surface, where a net 
uptake of O 2 from the atmosphere (20) and a net heating of the water occurs. The O 2 
flux/heating ratio in the Equatorial Pacific is thus opposite in sign to the cases considered above. 
This feature is a result of the upwelled waters having been exposed to the surface for a very brief 
period of weeks or days, which is insufficient to allow 0 2 to equilibrate with the atmosphere and 
for biological production to remove the nutrients. As the upwelled waters spread laterally away 
from the Equator, the net effect of warming and nutrient withdrawal leads to an overall O 2 
release to the atmosphere (20,28). Integrated over a wider latitude band, the net effect of 
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Equatorial upwelling on heat and 0 2 exchange is therefore more concordant with the main 
thermocline trend. 

The patterns noted above suggest the following generalizations. Although the changing 
efficiency of the biological pump dominates the 0 2 response of the ocean to warming and 
stratification, the 0 2 response is nevertheless strongly tied, over a range of space and time scales, 
to the net air-sea heat flux. For time scales of months to centuries, the 0 2 flux/heating ratios 
generally lie in the range of 2 to 10 nmol J^* . Larger ratios are found at higher latitudes, 
particularly in the Southern Hemisphere and for processes occurring over time scales of decades 
to centuries (e.g. thermocline ventilation) compared to time scales of months (seasonal 
exchanges). Some difference exist from ocean to ocean, with the North Atlantic having smaller 
0 2 response per unit heating or cooling than other high latitude regions. This analysis does not 
resolve the 0 2 response of the oceans to heating and cooling on thousand-year and longer time 
scales. 

ANTHROPOGENIC WARMING AND 0 2 OUTGASSING 

In Table 2 we use the 0 2 flux/heating relationships found above to formulate an estimate 
of the global air-sea 0 2 flux from 1990 to 2000. We consider three ocean regions: (1) North 
Atlantic (at all depths), (2) the deep Southern Ocean (> 1000m), and (3) the remaining oceans (at 
all depths). Warming in each region is assumed to produce 0 2 outgassing proportional to 
observed steady-state 0 2 */C' relationships in these regions. The approach is motivated by the 
observation that similar 0 2 /heat ratios are found for both transient warming on seasonal time 
scales and for steady-state warming/cooling over time scales from years to centuries. The 
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approach effectively adopts the "null hypothesis that the ocean s response to transient warming 
on decadal time scales is governed by similar ratios. 

Because compilations of ocean warming (1) have been completed only through year 
1998, and because these compilations do not resolve warming on a yearly basis below 300m 
depth, we rely on the model simulations (4,5) as a means of extrapolating the observed long-term 
( 1 950’s- 1 990's) warming in each region through the 1990-2000 period. The approach effectively 
accounts for the warming caused by anthropogenic radiative forcing but not natural variability. 
The latter we treat as a source of noise. We derive a central estimate of the global O 2 outgassing 
of0.29Pl0 14 moles 0 2 yr' 1 for 1990-2000, where 18% of the total 0 2 outgassing is attributable to 
warming in the North Atlantic, 20% in the deep Southern Ocean, and 61% in the remainder of 

the upper ocean. Our estimate corresponds to a global average 0 2 flux heating ratio of 5 nmol J 
1 

The uncertainty around our central estimate attributable to uncertainties in regional 
0 2 */& ratios and long-term warming rates is X0.13P10 14 moles 0 2 yr '. A much larger 
uncertainty must be allowed for decadal variability (15,29,30). One estimate of this can be 
derived by multiplying the global 0 2 /heat ratio of 5 nmol 0 2 J 1 by the decadal variability in 
global ocean heat storage, which we estimate from heat storage data (1), after removing the long- 
term trend, to be X5P10 22 J, which yields X2.5?10 14 mol 0 2 variability on a decadal basis. 
Here we adopt a slightly higher estimate of X4? 1 0 14 mol 0 2 , on the grounds that the heat 
storage data (1) may underestimate true variability due to spatial and temporal averaging, and 
given uncertainties in the appropriate 0 2 /heat ratio. Treating the decadal variability as a source 
of random noise, we derive an estimate of 0.29X0.4P 10 14 mol 0 2 yr for the total oceanic 0 2 
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outgassing from 1990 to 2000. 

From 1990 to 1998, the global upper ocean (< 300m) heat content increased at a rate of 
-0.5P10 22 J yr 1 (1), which is faster than the rate of -0.2P10 22 J yr" 1 which we would estimate 
from projecting long-term warming rates. If we assume that the additional warming of 
—0.3? 10 22 J yr'" 1 persisted though year 2000, and we add the difference to our outgassing 
estimate, then our central estimate increases from 0.29?10 14 mol yr 1 to 0.44P10 14 mol yr , 
assuming a scaling of 5 nmol J' 1 . To apply such a correction is premature, however, because 
the global 1990-1998 heating trend is heavily influenced by a large anomaly in the North 
Atlantic in 1998, which was possibly a transient associated with the 1997/98 El Nino event (1). 
In any case, the correction would be within our allowed uncertainties for decadal variability of 
X0.4P10 14 mol yr’ 1 . Once heat storage data are compiled through 2000, it may be possible to 
refine the decadal outgassing estimate and reduce the allowed uncertainties. 


CARBON BUDGETING 

What are the implications of oceanic 0 2 outgassing for carbon budgeting? A correction 
for oceanic 0 2 outgassing is needed to estimate land and ocean carbon sinks based on the global 
budgets of atmospheric 0 2 and C0 2 (3 1 ,32). These can be written, 

®C0 2 = F - O - B (2) 

®0 2 = -XfF + XbB + Z (3) 


where F is the C0 2 source from fossil-fuel burning, O and B are respectively the ocean and land 


12 


sinks, X F and X B are global-average 0 2 :C exchange ratios for fossil-fuel and land biota, and Z 
is the net source of 0 2 from the oceans. Equation (1) and (2) are solved to yield estimates of O 
and B, using (D0 2 and ®C0 2 derived respectively from observed changes in atmospheric 0 2 /N 2 
ratio and C0 2 mole fraction, F and X F derived from industrial records, and X B ^Dl.l (31,32). 

Previously it has been assumed that the ocean outgassing term Z is zero to within the 
uncertainties (31,32), or allowance has alternately been made for 0 2 outgassing based on the 
solubility effect alone (33). In Table 3, we correct the estimate of Manning (33), as cited in the 
recent IPCC (34) report, based on our above estimate of Z = 0.29X0.4 mol 0 2 /yr, where this 
estimate implicitly allows for the effect of anthropogenic warming on both solubility and 
stratification. The correction increases the oceanic sink by 0.18 Pg C yr' 1 and decreases the land 
sink by the same amount relative to the Manning (33) estimate. Although the change is small 
relative to other uncertainties, it nevertheless helps to reconcile the estimated oceanic sink with 
recent model estimates (35). 

A correction for oceanic 0 2 is also needed for carbon budgets based on direct 
measurement of changes in dissolved inorganic carbon (DIC) in the ocean. Detecting the input 
of anthropogenic C0 2 into the oceans is difficult due to large natural variability in DIC caused by 
ocean biology. This difficulty is commonly overcome by normalizing to a constant 0 2 
concentration to filter out the variability due to ocean biology (36,37,38,39,40). In effect, what 
is reported is not the change DIC, but rather the change in the quantity DIC + 0 2 /X, where X 
1.3 (35,36). By design, this approach neglects changes in DIC caused by variations in the 
efficiency of the biological pump. If the dissolved 0 2 inventory decreases globally due to 
stratification, the approach will underestimate oceanic C0 2 uptake by an amount given 
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approximately by (dI/dt)/X, where dl/dt is the change in oceanic 0 2 inventory. The correction 
is the same sign and a similar magnitude to that required based on atmospheric 0 2 and C0 2 
budgets. Taking the observed ocean warming of ~2 P 1 0 23 J between the middle 1950 s and 
1990's (1), and assuming the 0 2 inventory decreases by 5 nmol 0 2 J \ yields an upwards 
correction of ~9 Pg C over 1955-1995 period. In comparison, the total oceanic uptake from 

preindustrial times through 1990 is estimated to be 107X27 Pg C yr 1 (34). 

Recent estimates of change in oceanic inventories of DIC (38,39,40) have not allowed for 
changing 0 2 inventory, and therefore are presumably biased low, although more work is needed 
to establish reliable corrections in the individual ocean basins and globally. If the DIC data is 
additionally normalized based on regressions against £7 and alkalinity, then additional 
corrections may be needed for changing ocean heat content (independent of the effect on 0 2 ) and 
alkalinity inventory. 

DISCUSSION 

Our estimate of 0.29? 10 14 mol 0 2 yr' 1 for recent oceanic 0 2 outgassing caused by 
changing radiative forcing lies within the range estimated based on recent OGCMs (12,13,14,1 5). 
The relatively good agreement between these independent estimates strengthens the case that a 
long-term ocean outgassing of this magnitude may actually be occurring. Nevertheless, the 
similarity in these estimates undoubtedly results partly from similar assumptions. For example, 
the estimates commonly assume that Redfield P/0 2 ratios are constant with time and space, and 
that phosphate is the ultimate limiting nutrient. Also, while our estimate assumes that the 0 2 
response to transient warming is similar to the 0 2 response to steady-state warming and cooling, 
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the OGCMs similarly assume that the response of the ocean biology to transient warming can be 
understood based on parameters adjusted to duplicate steady-state behavior. It is unclear if these 
assumptions form a realistic basis for projection. 

Our estimate effectively assumes that O 2 */& relationships remain relatively constant 
during transient warming, or equivalently that preformed phosphate// 7 relationships remain 
constant during transient warming. For example, the O 2 flux/heating ratio we adopt for warming 
of the main thermocline would be obtained if the main effect of warming was to deepen the 
thermocline without altering the preformed phosphate/O relationship across the thermocline, or 
if the surface isotherms progress polewards during the transient at the same rate as the isolines of 
surface phosphate. The 0 2 flux/heating ratio we adopt for the deep Southern Ocean would be 
obtained if warming caused a reduction in deepwater formation around Antarctica, thus 
increasing the influence of other source waters on the chemistry of the deep ocean. 

Our estimate of the effect of Southern Ocean warming on 0 2 exchange can be compared 
to the study Broecker et al. (41), who used P0 4 \ 14 C, and chloro-flurocarbon (CFC) data to 
support the argument that the ventilation rate of deep waters around Antarctica slowed 
substantially in the 20 ,h century. According to their scenario, the PO 4 content of the deep ocean 
should now be decreasing with time, in which case a net 0 2 outgassing should also be occurring, 
as required by the links between PO 4 (i.e. 0 2 ) and the air-sea 0 2 flux. Taking their estimate of 
the difference in P0 4 * content of southern source waters relative to average deep waters and their 
estimate that the input rate has slowed by 1 0 ? 1 0 6 m 3 s^ 1 , yields a required outgassing rate of 
0 . 3 ? 10 14 mol 0 2 yr' 1 , which is five time larger than our estimate of 0.06P10 14 mol 0 2 yr ' 1 for the 
deep Southern Ocean (Table 2). We defend our smaller estimate on the following grounds: 
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First, the Broecker et al scenario implies that the deepwater heat content should be increasing by 
-0.4P10 22 J yr' 1 , based on the ~3“®C difference in temperature between southern surface waters 
and average deep waters. A warming rate of this magnitude is inconsistent with the observed 
temperature trends (1), unless the warming is mostly confined below 3000m, the maximum 
depth considered by Le vitus et al. (1). Second, recent work (42) suggests that a large 20 ,h 
century slow-down is not necessary to explain the CFC and PO 4 data. 

A reliable assessment of the global air-sea O 2 flux will ultimately require an approach 
based on direct observations rather than model studies. Over the next few decades, global 
systematic decreases in dissolved O 2 can be expected at the level of 0.4? 10 14 mol 02yr 1 or 
larger, which corresponds to a change of 0.7 Hmol O 2 kg^ 1 per decade, if spread uniformly over 
2000m. Resolving these changes against natural variability will require a high measurement 
density, but this may be feasible with appropriate sensor and platform development and with a 
concerted long-term observing program. Resolving these changes is needed for carbon 
budgeting, as a complement to lower-density DIC measurements and atmospheric O 2 /N 2 
measurements, and would help to assess the overall impact of climate change on the 
biogeochemistry and biodiversity of the oceans. Another source of uncertainty in carbon budgets 
involves change in oceanic organic carbon (18), which we have assumed is small, but is not well 
constrained. We suggest that a program to directly monitor oceanic inventories of O 2 and 
organic carbon, along with inorganic carbon and nutrients, should be given some priority in 
future ocean observing systems. 
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FIGURE CAPTIONS 

Figure 1. Ratio of the seasonal O 2 outgassing to the seasonal upper ocean warming by 10"® 
latitude bands. Dashed line indicates the ratios expected based on changes in O 2 solubility alone. 
Ratios from 20 = ®N to 20"®S are poorly defined due to weak seasonal forcing. The seasonal O 2 
outgassing here accounts for the anomalous flux after subtracting the mean annual outgassing. 

Figure 2. Scatter plot of the tracer O 2 * = 02 + 175 P0 4 versus potential temperature for data 
collected on the Geochemical Ocean Section Study (22, 23, 24). The data through the main 
thermocline of the Indian and Pacific oceans are seen to scatter around a line with a slope of 
approximately ^22Hmol kg 1 *® C^ 1 . Changes in O 2 are generally larger than expected from 
changes in the O 2 solubility (dashed line). Antarctic Surface Water from ref. 25. NADW: 


North Atlantic Deep Water. 
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Table 1. Significant recent changes in dissolved O 2 concentrations. 


Study 

Location 

Time 

Span 

Depth 
Range, m 

0 2 change 

Bmol kg' 

1 

Garcia et al. 
(10) 

Atlantic 24 section, 
averaged over full section 

1981- 

1992 

800-1900 

-3 to -7 

Pahlow and 
Riebesell (11) 

North Pacific, basin average 

1950's- 

1990’s 

unreported 

—5 

Keller et al. (9)‘ 

North Pacific, basin average 

1970's- 
1 990's 

300-500 

-1 to -2 




1000-1600 

+2 to +3 

Shaffer et al. 

(7) 

Pacific, 28"®S section, 
averaged from 52°®- 88°®W 

1967- 

1995 

800-1200 

-5 to -8 

Bindoff and 
McDougal (6) 

Indian Ocean, 32°®S section, 
averaged from 30‘®-117'®E 

1962- 

1997 

300-800 

-7 to -8 




2500-4000 

+3 

Matear et al. 
(13) 

Southern Ocean, 1 10°®- 
170"®E, 50°®-60"®S 

1965- 

1995 

>400 

-5 to -15 

Emerson et al. 
(8)’ 

North Pacific, 154“® W 
section, averaged from 22°®- 
44“®N 

1980- 

1997 

100-600 

-9 to -20 


‘Report change in O 2 -O 2 0 , rather than 0 2 concentration, where 0 2 ° is the solubility. 


Table 2. Estimated oceanic 0 2 outgassing 1990-2000 based on changing radiative forcing . 


Ocean Region 

Heating, 
10 22 Jyr' 1 

Heating, 

percent 

0 2 /heat ratio, 0 2 outgassing, 
nmol X 1 10 14 molyr' 

Outgassing, 

percent 

North Atlantic 

0.18X0.0 

8 

0.06X0.0 

n 

30% 

3X1.5 

0.05X0.04 

18% 

Deep Southern Ocean 

10% 

10X5.0 

0.06X0.04 

20% 

Other Oceans 

3 

0.36X0.1 

6 

60% 

5X2.5 

0.18X0.12 

61% 

Global Total 

0.6X0.18 

100% 

4.9 

0.29X0.13 

100% 
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'Heating rates based on observed warming in each region from 1955-1995 (1), scaled to yield 
total of 0.6? 10 22 J yr' 1 , the warming attributable to increasing radiative forcing from 1990-2000 
(4, 5). Error analysis assumes 45% uncertainty in individual heating rates and 50% uncertainty 
in 0 2 /heat ratios, with all errors uncorrelated. 
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Table 3. 02-based global carbon budget 1990-2000 (Pg C yr^ 1 ) . 



Manning (33) 

Revised budget 

Fossil-fuel burning 

6.33X0.4 

6.33X0.4 

Atmospheric C0 2 increase 

3.21X0.1 

3.21X0.1 

Total Ocean C0 2 sink 

1.68X0.5 

1.86X0.6 

Total Land C0 2 sink 

1.44X0.7 

1.26X0.8 


‘Budgets based on the decadal trend in atmospheric 0 2 /N 2 ratio at two northern hemisphere 
stations in the Scripps 0 2 /N 2 network.. Manning (33) assumes Z = 0.16X0.16?10 14 mol 0 2 yr' ; 
revised budget assumes 0.29X0.4? 10 14 mol 0 2 yr'‘(see text). Budgets differ in assumed ocean 
warming rate, and in assumed oceanic 0 2 outgassing per unit wanning. Both budgets allow for 
outgassing of N 2 (which impacts the observed 0 2 /N 2 ratio) at a rate based on N 2 solubility of 2.2 
nmol N 2 per J of warming. 







